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A&STRACT 


The  mathematical  models  of  thermoregulation  of  Stoiwijk  and  Hardy  and 
Montgomery  were  used  to  develop  a  model  suitable  for  the  simulation  of  human 
physiological  responses  to  cold-water  immersion.  Data  were  obtained  from 
experiments  where  thirteen  healthy  male  volunteers  were  totally  immersed 
under  resting  and  nude  conditions  for  1  h  in  water  temperatures  of  20  and  28°C. 
Mean  measured  rectal  temperature  (Tre)  fell  by  about  0.9  and  0.5° C  in  20  and 
2d°C  water  for  all  subjects,  yet  mean  measured  metabolic  rate  (M)  rose  by  aoout 
275  and  90  W  for  the  low  oody  fat  group  (n=7)  and  195  and  45  W  for  the  moderate 
body  fat  group  (n=6).  To  predict  the  observed  Tre  and  M  values,  the  present 
model  a)  included  thermal  inputs  for  shivering  from  the  skin  independent  of  their 
inclusion  with  the  central  temperature  to  account  for  the  observed  initial  rapid 
rise  in  M,  b)  determined  a  thermally  neutral  body  temperature  profile  such  that 
the  measured  and  predicted  initial  values  of  Tre  and  M  were  matcned,  c) 
confined  the  initial  shivering  to  the  trunk  region  to  avoid  an  overly  large 
predicted  initial  rate  of  rectal  cooling,  and  d)  calculated  the  steady-state 
convective  heat  loss  by  assuming  a  zero  heat  storage  m  the  skm  compartment  to 
circumvent  the  acute  sensitivity  to  the  small  skin-water  temperature  ^..'ference 
when  using  conventional  methods.  The  last  three  modifications  are  unique  to 
thermoregulatory  modelling.  A  UASIC  computer  listing  of  the  mocVi  and  a 
sample  simulation  are  provided. 


INTRODUCTION 

Nude  immersion  in  water  colder  than  the  deep  body  body  temperature 
represents  an  acute  exposure  to  cold  since  convective  heat  loss  is  many  times 
greater  than  in  air.  Temperature  gradients  becotne  large  and  physiological 
responses  are  dramatic.  These  factors  contribute  to  the  complexity  of 
mathematically  modelling  the  human  thermoregulatory  response.  Since  the 
inception  of  mathematical  models  of  human  thermoregulation  (see  reviews  by 
Hardy  (15)  and  Hwang  and  Konz  (IS)),  data  to  test  these  models  for  cold-water 
immersion  have  been  available  yet  have  been  only  recently  applied. 
Mathematical  models  of  thermoregulation  can  be  steady-state  or  dynamic. 
Steady-state  models  apply  where  a  heat  balance  exists,  and  therefore,  are 
limited  to  the  prediction  of  physiological  responses  that  do  not  change  with  time. 
Most  studies  of  cold-water  immersion  are,  however,  concerned  with  the  transient 
responses  upon  immersion.  Dynamic  models  can  be  applied  to  predict  these 
responses. 

Dynamic  models  of  thermoregulation  use  physical  representations  of  the 
human  body,  principles  of  heat  conduction,  and  control  theory  to  simulate 
physiological  responses  to  a  change  in  the  environment.  Such  models  provide  a 
useful  theoretical  device  to  evaluate  and  interpret  experimental  data,  and 
potentially  can  be  applied  to  a  wide  range  of  subject  classifications  and 
environmental  conditions.  The  dynamic  models  assessed  by  Hardy  in  1972  (15) 
were  found  inadequate  (or  predicting  human  responses  to  cold  environments. 
Among  these  was  the  Stofwtjk-Hardy  model  (25)  o.igmaUy  developed  to  predict 
the  physiological  responses  of  nude  man  in  an  air  environment.  In  1976,  Gordon 
et  at.  (13)  extended  the  concepts  of  Stolwijk  and  Hardy  (25)  and  'insslor  (29)  to 
model  the  physiological  responses  to  a  transient  cold  air  exposure.  In  1QS4, 
Wissler  (30}  evaluated  this  model's  predictive  capability  for  cold-water 
immersion  and  found  it  unsuitable.  Other  models  that  Wisder  evaluated  included 


his  own  and  Stolwijk-Hardy,  yet,  the  agreement  between  measured  and  predicted 
values  of  temperature,  metabolic  rate,  and  net  sensible  heat  loss  was  found  to  be 
less  than  satisfactory.  Although  a  more  recent  application  of  the  Wissler  model 
to  other  cold  water  immersion  studies  has  provided  improved  predictions  (3  a 
difficulty  with  this  model  is  its  inability  to  match  measured  and  predicted  initial 
core  temperatures  and  metabolic  rates.  Strong  and  Goldman  (28)  developed  a 
linearized  model  for  predicting  skin  and  rectal  temperatures  specifically  for  cold 
water  immersion.  However,  their  model  did  not  include  conductive  heat 
exchange  between  tissue  and  blood  which  is  known  to  critically  affect  the  heat 
storage  of  the  body,  nor  did  they  model  for  a  muscle  compartment  where 
shivering  activity  increases  the  body’s  metabolic  rate  (the  model  requires  the 
experimentally  determined  metabolic  rate  as  an  input). 

An  alternative  model  is  a  version  of  the  Stolwijk-Hardy  model  developed  by 
Montgomery  (£1)  also  tor  coid  water  immersion,  but  not  tested  by  Wissler  (30). 
In  our  study  using  data  of  resting  nude  subjects  totally  immersed  in  cold  water, 
we  found  the  Montgomery  model  not  wholly  satisfactory  for  predicting  transient 
changes  in  rectai  temperature  and  metabolic  rate,  yet  with  certain  modifications 
of  the  model,  good  agreement  was  obtained.  These  modifications  include  adding 
a  shivering  component  responsive  to  skin  temperature  only,  matching  the 
measured  and  predicted  initial  core  temperatures  and  metabolic  rates,  Confining 
the  initial  shivering  to  the  trunk  region,  and  determining  the  steady-state 
convective  heat  toss  to  the  water  through  use  of  the  lieat  storage  eauation. 
Various  forms  of  the  first  modification  have  appeared  elsewhere  m  other  models 
whereas  the  lost  three  modifications  are  unique  to  thermoregulatory  modelling 
and  may  be  generally  applicable. 

This  report  describes  these  modifications  and  presents  a  comparison 
between  the  measured  and  predicted  thermoregulatory  responses  for  whole  body 
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immersion  in  cold  water.  Data  from  experiments  of  nude  whole  body  immersion 
were  used  since  skin  temperatures  quickly  reach  steady-state  values  and  the  heat 
losses  to  the  water  are  limited  to  convective  heat  transfer  which,  as  will  be 
seen,  can  easily  be  determined  once  a  steady-state  skin  temperature  is  reached. 
A  computer  listing  of  the  model  and  a  sample  simulation  are  given  in 

APPENDIX  J 

METHODS 

Data 

Data  presently  used  were  available  from  a  series  of  whole  body  water 
immersion  studies  (5,12).  Thirteen  healthy  male  volunteers  were  totally 
immersed  under  resting  and  nude  conditions  for  1  h  in  water  temperatures  of  20 
and  2S°C.  Since  a  subject's  thermal  and  metabolic  response  depends  largely  on 
his  body  composition  (S,17,19,2Q,24,2S),  this  study  classified  the  subjects  into 
two  groups,  those  of  low  body  fat  (LBF)  and  those  of  moderate  body  fat  (MBF). 

Mean  (♦  SO)  anthropometric  values  of  the  LBF  group  (n^7)  were: 
height  =  174.9  (4.5)  cm,  weight  -  69.G  (7.5)  kg,  skinfold  =  5»$3  (Q.7y)  mm,  body 
fat  s  *7.67  (1.57)%,  and  surface  area  :  l.X)  (0.11)  tn^.  Values  for  the  MDF  group 
(n«b)  were:  height  =  175.7  (6.4)  cm,  weight  :  79.2  (13.1)  kg,  skinfold  *  H.S2 
(4.26)  mm,  body  fat  2  17.62  (4.11)%,  and  surface  area  =  1.94  (0.17)  m*. 

Model 

The  model  used  m  the  present  study  is  based  largely  on  the  Montgomery 
version  (21)  of  the  StoJwijk -'Hardy  model  (25).  The  human  body  is  treated  as  a 
passive  heat  transfer  system  and  is  divided  into  six  d’stmct  segments,  the  head 
modelled  as  a  sphere  and  the  trunk,  arms,  h^tvS,  legs,  and  feet  modelled  as 
cylinders.  The  model  is  shown  schematically  m  Fig.  i.  Heat  flows  radially  in  the 
model  segments  and  heat  transfer  between  segments  is  through  conduction  via 


the  central  blood.  Each  segment  is  composed  of  four  concentric  annular 
compartments,  the  core,  muscle,  fat,  and  skin,  as  proposed  by  Stolwijk  and  Hardy 
(25).  In  addition,  the  central  blood  is  a  single  compartment  located  within  the 
trunk  segment.  The  Montgomery  version  of  expanding  the  number  of  core  an>; 
muscle  compartments  by  four  each  is  not  used.  Instead,  the  relative  weight 
distribution,  thermal-capacitance  values,  basal  metabolic  rates,  and  basal  blood- 
flow  rates  of  all  compartments  proposed  by  Montgomery  are  used.  The  reader  is 
referred  to  Ref.  21  for  these  values. 

The  thermoregulatory  controlling  system  integrates  the  thermoreceptf' 
output  signals  of  certain  compartments  and  determines  the  response  th.  igh 
efferent  commands.  For  example,  cold  signals  may  induce  shivering.  The 
thermoreceptor  output  signal  of  each  compartment  is  determined  by  tne 
difference  between  the  compartment's  current  temperature  and  its  set-point 
value.  Set-point  values  are  established  before  immersion  and  remain  constant 
throughout  the  immersion.  The  efferent  commands  involve  sweating,  vasomotor 
response,  and  shivering.  Unless  otherwise  indicated,  the  thermoreceptor  and 
efferent  output  simulations  follow  tire  method  of  Montgomery  (21). 

Thermal  conductances  between  compartments  were  determined  according 
to  the  method  outlined  by  btolwijk  and  Hardy  (2».  Thermal  resistances  for 
spherical  and  cylindrical  geometries  were  obtained  from  Sekms  and  Emery  \21). 
Thermal  conductivity  values  for  the  core  and  muscle  compartments  were  taken 
from  Stoiwifk  (2b),  and  those  for  the  fat  and  skm  compartments  were  taken  from 
Sekms  and  Emery  (2J). 

Since  the  subjects  were  totally  immersed  m  water,  both  radiative  and 
evaporative  heat  transfer  from  the  body  were  considered  negligible.  Total 
respiratory  heat  loss  was  determined  by  the  combined  respired  evaporative  and 
respired  convective  heat  losses  of  the  trunk  core  ii!)  and  by  the  basal 


M 


respiratory  heat  loss  of  head  core  (26).  The  subjects  breathed  through  a  snorkel 
and  therefore  the  respiratory  heat  loss  was  determined  by  assuming  that  the  air 
breathed  was  fully  saturated  and  at  a  temperature  equal  to  the  water 
temperature. 

Initial  and  Set-Point  Temperatures 

Initial  conditions  assume  thermal  neutrality.  By  simulating  an  exposure  to 
an  arbitrary  environment  in  the  zone  of  thermal  neutrality,  the  original  Stolwijk- 
Hardy  model  will  generate  steady-state  temperatures  for  all  compartments 
including  the  central  blood  (26,27).  The  resultant  initial  temperatures  are  thus 
assigned  as  the  set-point  temperatures  for  thermoregulation.  One  drawback  with 
this  method  is  that  the  model's  temperature  profile,  which  is  based  on 
"standard"  man,  does  not  necessarily  match  the  subject's  profile  in  his  pre¬ 
immersion  state,  and  therefore,  thermoregulation  may  be  arbitrarily  imposed. 
Furthermore,  initial  offsets  between  measured  and  predicted  core  temperatures 
may  affect  the  level  of  agreement  during  the  subsequent  immersion  phase. 

Ideally,  the  initial  temperature  profile  of  the  model  should  march  the 
subject's.  At  present,  it  is  not  possible  to  measure  the  subject’s  temperature 
profile,  so  certain  assumptions  must  be  made.  First,  in  accordance  with  Stolwijk 
(26),  it  was  assumed  that  the  subject  was  thermally  neutral  in  his  pre-immersion 
state  (this  is  reasonable  considering  that  the  subjects  in  our  study  were  resting  in 
an  air  environment  within  the  /.one  of  thermal  neutrality  before  immersion). 
Second,  it  was  assumed  that  the  subject's  measured  pre-immersion  metabolic 
rate  represented  his  basal  value  (BMK)  and-  that  the  model's  trunk  core 


temperature  represented  his  rectal  temperature  (Tre). 


A  thermal  neutral  temperature  profile  can  thus  be  determined  for  any 
subject  by  setting  the  heat  storage  of  each  model  compartment  equal  to  zero  and 
solving  the  resulting  linear  equations  using  matrix  methods.  By  specifying  values 
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of  the  trunk  core  (rectal)  and  central  blood  temperatures,  an  iterative  solution  is 
sought  such  that  the  heat  storage  of  the  central  blood  is  also  zero.  The 
temperature  distribution  obtained  for  "standaid"  man  by  this  procedure  agrees 
exactly  to  that  obtained  using  the  convention  procedure,  but  for  conditions  other 
than  "standard",  this  procedure  has  the  advantage  of  matching  model  and  subject 
values  of  Tre  and  BMR.  The  predicted  neutral  temperature  profile  obtained 
using  this  procedure  for  the  MBF  group  is  shown  in  Fig.  2. 

Convective  Heat  Loss 

A  major  theoretical  obstacle  for  any  thermoregulatory  model  is  the 
determination  of  convective  heat  loss,  especially  in  water  immersion  where  heat 
transfer  is  many  times  greater  than  in  air  (4,24).  Heat  transfer  is  sensitive  to 
the  skin-water  temperature  difference,  especially  as  the  skin  temperature  (T^) 
nears  the  temperature  of  the  water.  Because  of  this  sensitivity,  small  changes  in 
T^  can  cause  large  changes  in  the  predicted  convective  heat  transfer  (C),  as 
demonstrated  in  the  DISCUSSION.  These  large  changes  in  C  critically  affect  the 
heat  st-.*  age  of  the  skin  compartment  (S^),  and  consequently  the  heat  loss  of  the 
oody.  This  problem  is  exacerbated  by  the  assumptions  of  body  shape  and  water 
motion  that  determine  the  heat  transfer  coefficient  (see  APPENDIX  A). 

It  has  been  shown  experimentally  that  the  mean  weighted  skin  temperature 
(T^)  of  nude  subjects  falls  exponentially  during  immersion  in  cold  water  (28), 
and  that  the  asymptotic  limit  is  a  temperature  slightly  higher  than  the  water 
temperature,  although  the  skin-water  temperature  difference  increases  with 
lowered  water  temperature  (20,22,28).  These  experimental  observations  can  be 
coupled  with  the  theoretical  determination  of  the  convective  heat  transfer 
coefficient  to  arrive  at  a  model  prediction  of  convective  heat  loss  that  avoids 
the  uncertainties  discussed  above. 
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First,  the  rate  of  change  of  skin  temperature  (i^)  must  equal  Sg^/C^ 
where  Cg^  is  the  heat  capacity  of  the  skin.  Second,  Tg^  should  be  proportional  to 
the  difference  between  the  skin  compartment's  steady-state  temperature  (T^  ) 
and  its  current  temperature  to  approximate  the  exponential  fall  in  Tg^  and  to 
allow  Tgk  to  approach  T,^  asymptotically.  Through  numerical  integration,  the 
incremental  change  in  Tgg  can  thus  be  approximated  by: 

AT"sk  =  (Tskss  -  Tgj^)  [1  -  exp  (-Sg^*  At/Cg^  (Tskss  -  Ts|<o))j,  (1) 

where  Tg^  is  the  skin  temperature  before  the  incremental  change  and  At  is  the 
time  increment  chosen  sufficiently  small  so  that  the  above  constraints  are 
satisfied.. 

In  the  present  study,  T^  was  assigned  the  experimentally  measured 
value,  yet  an  arbitrary  value  close  to  the  temperature  of  the  water  could  have 
been  assigned  without  incurring  a  large  error  in  determining  the  convective  heat 
loss  (see  DISCUSSION).  The  heat  storage  of  the  skin  was  determined  through  a 
thermal  balance  of  the  skin  compartment: 

^sk  ~  Msk  ~  ~  '  ^skol  +  Kfsk>  (2) 

where  Mg^  is  the  metabolic  rate  of  the  skin,  Kg^bj  is  the  conductive  heat 
transfer  rate  from  the  skin  to  the  blood,  and  Kfgj<  is  the  conductive  heat  transfer 
rate  from  the  fat  to  the  skin.  The  convective  heat  transfer  was  determined 
through  fluid  dynamic  considerations  (see  Eq.  6  and  APPENDIX  A).  This 
calculation  was  carried  out  by  assuming  a  water  velocity  of  0.005  m/s  which 
represents  the  motion  produced  in  "still"  water  by  respiration  and  mild  shivering 
(32).  Although  shivering  intensity  can  be  expected  to  increase  with  increased 
immersion  time,  a  steady-state  skin  temperature  was  attained  weli  before  the 
water  motion  was  seriously  underestimated. 

Once  the  skin  temperature  was  close  to  its  assigned  steady-state  value 
(assumed  by  the  model  when  the  difference  between  Tgg  and  Tskss  'vas  Iess  than 
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O.QO.WC),  no  further  change  in  skin  temperature  occurred.  The  convective  heat 
transfer  from  skin  to  water  was  then  determined  assuming  zero  heat  storage  of 
the  skin  compartment  (i.e.  setting  =  0  in  Eq.  2). 

Efferent  Shivering  Command 

Central  to  any  thermoregulatory  model  for  cold  exposure  is  the  efferent 
command  for  shivering.  Montgomery  (21)  specified  the  shivering  command  as  a 
product  of  a  control  coefficient,  a  central  (head  core)  thermoreceptor  output 
signal  and  the  appropriate  skin  (peripheral)  thermoreceptor  output  signal.  As 
will  be  seen,  such  an  expression  is  incapable  of  predicting  the  initial  rapid  rise  in 
metabolic  rate  that  has  been  repeatedly  observed  for  cold-water  immersion 
(1,10,16).  There  is  sufficient  evidence  to  support  the  view  that  to  some  extent, 
shivering  is  independently  controlled  by  skin  thermoreceptors  (1,3,6,8,14,28).  In 
fact,  the  original  Stolwijk-Hardy  model  allowed  for  this.  Since  the  initial  rapid 
incre-je  in  metabolic  rate  correlates  well  with  the  observed  initial  rapid 
decrease  in  skin  :  emperature,  the  controller  equation  for  shivering  in  the  present 
model  included  a  shivermg  component  responsive  to  skin  temperature  only. 

An  increase  in  the  metabolic  rate  due  to  shivering  entails  a  corresponding 
increase  in  muscular  bl-.od  ilow  which  can  indirectly  affect  the  core 
temperature.  If  the  arm  and  leg  muscle  temperatures  are  lower  than  that  of  the 
central  blood  as  indicated  in  Fig  2,  then  any  sudden  increase  in  blood  flow  to 
these  muscles  will  lower  the  central  bbod  temperature.  Given  that  the  observed 
metabolic  rate  initially  rises  rapidly,  a  model  prediction  of  a  cor-esponding 
increased  biood  flow  to  the  limb  muscles  would  indirectly  cause  an  initial  fall  in 
trunk  core  temperature  (through  conductive  heat  exchange  with  the  central 
blood)  much  mot*'  rapidly  than  observed.  To  avoid  this,  the  present  model 
confined  initial  shivering  to  the  trunk  (since  its  temperature  was  close  to  that  of 
the  central  blood)  and  delayed  ♦he  onset  of  shivering  of  the  limb  muscles 
exponentially. 
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Counter-Current  Heat  Exchange 


To  conserve  body  heat,  counter-current  heat  exchange  of  the  limbs  (i.e. 
arms,  hands,  legs,  and  feet)  may  occur.  A  simple  yet  effective  means  of 
modelling  this  is  to  assume  an  effective  temperature  of  the  blood  in  the  limb 
which  is  used  to  determine  the  conductive  heat  exchange  in  that  compartment 
temperature  (in  which  case  the  counter-current  heat  exchange  is  100% 
effective).  The  expression  used  to  determine  this  value  is 

TBL(i)  =  T(i)  +  (Tbl  -  T(i)>exp  (-X*  COLDS)  (3) 

where  TBL(i)  is  the  effective  blood  temperature  in  the  ith  compartment,  T(i)  is 
the  compartment  temperature,  T51  is  the  central  blood  temperature,  X  is  a 
proportional  control  coefficient,  and  COLDS  is  the  weighted  skin  thermoreceptor 
output  signal.  The  dependence  on  COLDS  allows  for  an  increase  in  the  counter- 
current  heat  exchange  with  increasing  severity  of  exposure. 

Simulation  Procedure 

The  anthropometric  characteristics  of  the  model  subject  assumed  the 
average  values  for  the  group  it  was  simulating.  The  neutral  (and  set-point) 
temperature  profile  was  determined  separately  for  each  body  fat  group  and 
exposure  based  on  the  group's  mean  measured  pre-immersion  Tre  and  M.  Values 
of  air  temperature  (Tair_neutrd|)  and  central  blood  temperature  determined  for  a 
condition  of  thermal  neutrality  are  listed  m  Table  l.  During  the  immersion,  the 
compartments'  heat  storage  were  determined  using  the  finite  difference 
procedure  outlined  by  Stolwijk  (26).  The  incremental  change  m  temperature  of 
any  compartment  could  not  exceed  0.1  ^C. 


RESULTS 

Figures  3  to  6  illustrate  the  measured  (USE)  and  predicted  values  of  the 
rectal  (modelled  as  the  trunk  core)  temperature  and  metabolic  rate.  To  obtain 
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these  predicted  values,  the  following  controller  expression  for  shivering  was 
used*. 

CHILL  =  AD*l5-COLl)(l)-COLDS  ♦  65*(COLDS/PBF)L5),  (4) 

where  CHILL  is  the  metabolic  response  (W)  to  the  cold  stress,  Aq  is  the  subject's 
surface  area  (m^),  COLD(i)  is  the  head  core  thermoreceptor  output  (equal  to  the 
difference  between  the  current  temperature  of  the  head  core  and  its  set-point 
value  only  when  the  head  core  temperature  is  less  than  its  set-point  value, 
otherwise  the  output  value  is  zero),  and  PBF  ts  the  subject's  percent  body  fat. 
For  the  MBF  group  immersed  in  28°C  water,  a  value  of  2  instead  of  5  was  used 
for  the  proportional  control  coefficient  of  the  first  term.  In  all  cases,  the  value 
of  X  was  zero  indicating  that  it  was  not  necessary  to  use  the  counter-current 
heat  exchange  mechanism,  althougn  it  remains  in  the  model  as  an  option.  In  the 
Montgomery  model  (21),  only  the  first  term  of  Lq.  4  is  present,  and  the  product 
of  Aq  times  the  control  coefficient  was  assigned  a  value  of  24.4  W  (21  kcal/h). 

To  avo»d  the  excessive  initial  decrease  in  trunk  core  temperature  discussed 
earlier,  initial  shivering  was  confined  to  the  trunk  and  a  portion  of  this  shivering 
was  shifted  to  the  arm  and  leg  muscles  exponentially  according  to: 

CHILM  (trunk)  Q.S)  ♦  0. 12  exp  (-0.W/PBF) 

CHILM  (arm)  ■  0.05  (1-exp  (O.H/PBf *))  (3) 

CHILM  (leg)  h  0.07  (t-exp  (-0.W/PBP)), 

whore  CHILM  is  the  weighing  factor  of  the  corresponding  muscle's  contribution 
to  the  overall  shivering,  t  is  the  elapsed  time  (mm)  since  immersion,  and  the 
control  coefficients  Q.$3,  0.03  and  0.07  were  taken  from  Stolwijk  and  Hardy  (27). 
Note  that  as  t  increases,  the  CHILM  factors  revert  to  the  values  given  by 
Stolwijk  and  Hardy  (27)  which  were  also  used  by  Montgomery  (21).  The 
attenuation  by  the  group's  PBF  jf  both  the  shivering  command  in  Lq.  4  and  the 
exponent  in  Lq.  ^  was  a  necessary  modelling  construct  to  obtain  the  results 
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shown  in  Figs  3  to  6.  Also  shown  in  thes-i  figures  are  the  predicted  values  using 
the  present  model  but  without  the  independent  shivering  command  from  the  skin 
and  without  the  delayed  onset  of  shivering  of  the  limb  muscles. 

The  predicted  temperature  profile  for  the  MBF  group  after  1  h  of 
immersion  in  20°C  water  is  illustrated  in  Fig.  2.  Every  compartment  except  the 
trunk  muscle  shows  a  decrease  in  temperature;  the  increase  in  temperature  in 
the  trunk  muscle  is  slight,  from  37.10  to  37.33°C.  Decreases  in  temperature  in 
the  other  compartments  range  from  small  changes  in  the  core  and  muscle  of 
active  compartments  to  large  changes  in  the  inactive  compartments  and  the  fat 
and  skin  compartments  of  all  segments.  Figure  2  is  representative  of  the  model 
prediction  (in  a  qualitative  sense)  of  the  L13F  group  and  for  immersion  in  28°C 
water  of  both  groups. 

Figure  7  shows  the  model  prediction  of  mean  body  temperature  (T|j),  trunk 
core  temperature,  mean  skin  temperature,  metabolic  rate,  and  convective  heat 
loss  for  the  MISF  group  immersed  in  20°C  water.  The  mean  body  temperature 
was  determined  by  weighting  each  compartment's  temperature  according  to  its 
heat  capacity  (26),  The  mean  skin  temperature  (T^)  was  similarily  determined 
from  all  skin  compartments.  The  overall  convective  heat  loss  was  determined  by 
summing  the  convective  heat  loss  of  each  segment. 


DISCUSSION 


Efferent  Shivering  Command 

To  obtain  agreement  with  the  measured  metabolic  and  thermal  response  to 
cold-water  immersion,  an  efferent  shivering  command  based,  in  part, 
independently  on  the  skm  temperature,  and  a  delayed  onset  of  limb  shivering  was 


required.  The  possibility  of  an  independent  skin  temperature  effect  on  shivering 
was  not  excluded  in  the  Stolwijk-Hardy  and  .Montgomery  models  (although  it  was 


not  used  by  Montgomery  (21)),  and  as  pointed  out  by  Cabanac  (6),  the  debate  over 
additive  versus  multiplicative  combinations  of  thermoreceptor  output  signals  has 
not  been  resolved.  In  the  present  model,  the  shivering  command  from  the  skin 
appears  to  be  dependent  on  the  skin  thermoreceptor  output  signal  raised  to  the 
power  1.5.  Furthermore,  this  signal  is  attenuated  by  the  subject's  percent  body 
fat,  also  raised  to  the  same  power.  Differences  in  the  shivering  response  to  the 
same  core  and  skin  temperatures  between  low  and  moderate  body  fat  groups  has 
been  previously  reported  for  the  data  used  in  this  study  (2S). 

The  independent  efferent  command  from  the  skin  was  necessary  to  predict 
the  observed  initial  rapid  rise  in  metabolic  rate.  This  is  demonstrated  by  the 
dashed  lines  in  Figs.  3  through  6  where  only  the  first  term  of  Eq.  4  was  used  and 
its  coefficient  was  adjusted  to  correspond  to  the  value  used  by  Montgomery  (21). 
One  reason  that  the  rapid  initial  rise  in  metabolic  rate  cannot  be  predicted  with 
the  shivering  command  based  on  the  product  CCLlMD  •  CQLDb  (see  Eq.  4)  alone 
is  that  the  head  core  temperature  is  very  slow  to  change  initially  (1,3),  and 
therefore,  despite  the  rapid  initial  change  in  skm  temperature,  the  product  of 
cold  signals  from  the  head  core  and  skin  has  a  depressed  value  in  the  initial  stage 
of  immersion.  In  fact,  the  central  temperature  may  initially  increase  (10)  in 
which  case  the  product  has  a  /.eru  value. 

An  alternative  shivering  command  could  liave  been  based  on  the  time 
derivative  of  the  skm  temperature  (30).  Such  a  command  would  produce  a 
transient  increase  in  shivering  intensity.  Considering  the  rapidly  falling  Tjg  upon 
immersion,  this  transient  would  decay  well  before  any  appreciable  decrease  in 
Trc.  Instead,  our  data  indicated  that  initial  values  of  M  peaked  between  6  and  IS 
minutes  after  immersion,  much  longer  than  the  few  minutes  it  took  for  a  steady- 
state  skin  temperature  to  be  reached.  Because  of  the  high  variability  i r. 
individual  responses,  no  attempts  were  made  to  model  this  behavior. 


■p  delayed  onset  of  limb  shivering  was  necessary  to  avoid  a  model 
prediction  of  a  large  initial  decrease  in  trunk  core  temperature.  Such  a  decrease 
would  stem  from  increased  blood  flow  from  the  cooler  muscles  of  the  arms  and 
legs  thereby  lowering  the  central  blood  temperature  which  in  turn  would  lower 
the  trunk  core  temperature  (27).  The  exponential  factor  governing  the  delayed 
onset  (see  Eq.  5)  suggests  that  limb  shivering  of  the  LBF  group  began  sooner 
than  that  of  the  MBF  group.  At  present,  direct  experimental  evidence  to  test 
this  dependence  on  body  fat  is  lacking. 

Set-Point  Temperatures 

The  thermal  neutral  temperature  profile  and  hence  the  set-point  values  for 
thermoregulation  were  determined  according  to  the  pre-immersion  data  of  the 
subjects  and  not  on  the  expected  values  for  the  standard  man  as  used  in  the 
Stolwijk-Hardy  model  (26,27).  Th"  possibility  of  adjustable  set-point 
temperatures,  which  our  method  inherently  assumes,  has  been  reported 
previously  (14).  The  advantage  that  the  present  method  provides  over  the 
Stolwijk  method  is  to  assure  that  the  model  subject  is  thermally  neutral  at  the 
outset  of  an  exposure  and  that  the  measured  and  predicted  initial  values  of  core 
temperature  and  metabolic  rate  are  matched.  This  procedure  is  not  limited  to 
cold  water  immersion  and  may  be  potentially  useful  for  all  environmental 
conditions. 

Tissue  Conductance  and  Heat  Transfer  Coefficient 

An  important  test  and  useful  application  of  the  present  model  is  its 
prediction  of  average  tissue  conductance,  k,  and  the  convective  heat  transfer 
coefficient,  H^.  These  values  can  be  calculated  from  the  model  predictions  as 
(4): 

k  =  C/(Tre  -  f^).  (6) 

and 

»v  *  c/tr*  -  rw),  (7) 
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where  Tre  is  represented  by  the  trunk  core  temperature  and  Tw  is  the  water 
temperature.  Table  2  lists  these  values  for  both  body  fat  groups  and  exposures 
after  1  h  of  immersion.  Steady-state  conditions  can  be  assumed  at  this  time  (2), 
as  demonstrated  in  Fig.  7.  The  values  of  average  tissue  conductance  shown  in 
Table  2  are  in  good  agreement  with  other  reported  values  (4,10,20,24).  In  fact, 
the  predicted  increase  of  average  tissue  conductance  with  lowered  water 
temperature  is  consistent  with  the  decreasing  insulative  value  of  increasing!) 
active  muscle  (19).  Such  a  decrease  was  noted  by  both  Craig  and  Drovak  (10)  and 
McArdle  et  al.  (20)  where  Tw  was  lowered  from  23  to  24°C.  Further  support  of 
the  model  stems  from  its  prediction  of  higher  average  tissue  conductance  for  the 
LBF  group  compared  to  the  M15F  group. 

The  model-predicted  values  of  the  convective  heat  transfer  coefficient 
(see  Table  2)  are  in  agreement  with  the  values  reported  by  Witherspoon  et  al. 
(32),  Nadel  et  al.  (22),  and  Strong  et  al.  (23)t  but  are  much  higher  than  those 
reported  by  tioutelier  et  al.  (2).  The  potential  for  such  a  disparity  has  already 
been  noted  by  lioutelier  et  al.  (2)  and  reasons  given  stem  from  differences  in  ttie 
measurement  and  theoretical  determination  of  convective  heat  loss.  It  should  be 
noted  that  the  heat  transfer  coefficient  is  highly  sensitive  to  the  skin-water 
temperature  difference.  For  instance,  complete  agreement  between  the  he 
values  for  the  20  and  2$°C  exposures  of  either  body  fat  group  can  be  obtained  by 
increasing  the  steady-state  skin  temperature  by  less  than  0.1  °C  for  the  exposure 
to  2S°C  water. 

The  procedure  by  which  Ute  present  model  determined  the  convective  heat 
toss  to  the  water  avoided  this  sensitivity  once  steady-state  of  the  skin 
temperature  was  reached.  Recall  that  the  experimentally -measured  values  of 
the  steady-state  skin  temperature  were  used,  however,  the  choice  of  T ^  could 
have  been  made  arbitrarily  without  significantly  offset  mg  the  final  result  since 


the  convective  heat  loss  during  steady-state  of  the  skin  temperature  is  largely 
determined  by  the  conductive  heat  transfer  from  the  fat  to  the  skin.  This  heat 
transfer  is  only  slightly  affected  by  small  changes  in  the  steady-state  value  of 
the  skin  temperature.  For  example,  if  Tgp^  was  raised  from  21.0  to  21.5°C  for 
the  MBF  group  immersed  in  20°C  water,  the  predicted  convective  heat  loss  to 
the  water  would  change  by  less  than  4%  from  166.0  to  160.2  W/m^.  Note, 
however,  that  changing  the  skin  temperature  from  21.0  to  2i.5°C  would  decrease 
the  convective  heat  transfer  coefficient  (See  Eq.  7)  by  67%  which  further 
demonstrates  the  potential  disparity  among  reported  values  of  hc  as  pointed  out 
by  Bouteiier  et  al.  (2). 

Mean  Body  Temperature 

The  present  model  can  provide  insight  into  the  thermal  response  of  the 
whole  body  from  its  prediction  of  the  rate  of  change  of  mean  body  temperature. 
Assuming  that  the  mean  body  temperature  can  be  approximated  by  (9): 

Tb  =xTsk  »(l -x)Tre,  (8) 

where  x  varies  depending  on  the  environment,  then  the  rate  of  change  of  mean 
body  temperature  (Tb)  should  be  less  than  the  rate  of  change  of  rectal 
temperature  (T^e)  after  the  skin  temperature  has  reached  its  steady-state  value. 
Since  the  fecial  temperature  c  *n  represented  by  the  trunk  core  temperature, 

this  prediction  holds  true,  as  car.  be  seen  from  dv*  estimr'ed  (based  on  the  slope 

• 

of  temperature  against  time)  values  of  Tb  and  listed  in  Table 

To  check  on  the  internal  consistency  of  the  model,  the  rate  of  change  of 
mean  body  temperature  can  alternatively  be  determined  through  the  thermal 
balance  equation  for  a  body  totally  immersed  m  water  by  (21); 

Tb  (M-c:-HR)/Cb,  (9) 

where  HR  is  the  rate  of  total  respiratory  heat  loss  (from  head  and  trunk  core 
compartments)  and  Cp  is  the  heat  capacity  of  the  whole  body.  The  values  of  Tb 

1> 


determined  by  Eq.  9  and  shown  in  Table  2  are  in  close  agreement  with  those 
estimated  from  the  slope  of  temperature  change  with  time.  This  confirms  that 
the  model  is  self-consistent  with  the  prediction  of  changing  body  temperature 
during  cold-water  immersion. 

Conclusion 

Nude  immersion  in  cold  water  summons  dramatic  physiological  responses 
not  fully  considered  in  the  early  development  of  thermoregulatory  models.  To 
model  this  response  mathematically  requires  refinement  of  certain  mechanisms 
that  are  otherwise  adequate  for  less  acute  exposures.  This  was  the  rationale  for 
the  modifications  of  the  Stolwijk-Hardy  and  Montgomery  models  from  which  the 
present  model  evolved,  Although  these  modifications  were  derived  explicitly  for 
cold  water  immersion,  they  may  be  generally  applicable  to  other  conditions. 

The  inclusion  of  an  independent  shivering  command  from  the  skin  was 
optional  in  tlie  Stolwijk-Hardy  and  Montgomery  models  although  it  was  not  used 
by  Montgomery  (21).  Without  this  independent  shivering  command,  it  is  not 
possible  to  predict  the  initial  rapid  increase  m  metabolic  rate  for  nude 
immersion  m  cold  water  using  onl>  the  product  of  signals  from  the  head  core  and 
the  skin. 

the  remaining  modifications  are  unique  to  thermoregulatory  modelling. 
Die  prevent  method  of  determining  a  thermally  neutral  temperature  profile 
allows  matching  of  the  predicted  and  measured  initial  core  temperatures  and 
metabolic  rates.  Die  delayed  onset  of  hmb  shivering  eases  the  transition  of 
increased  blood  flow  to  tin*  limb  muse les  thereby  avoiding  too  rapid  a  decrease  in 
central  biood  ttmperatu? e.  The  use  of  the  heat  storage  equation  of  the  skin 
compartment  to  predict  the  convective  heat  loss  during  steady-state  of  the  skin 
temperature  circumvents  the  high  sensitivity  to  the  skin-water  tempe'atute 
difference  witen  using  conventional  methods. 
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TaDle  1:  Measured*  and  model  values  for  thermal  neutrality  in  air  and  for 


response  to  cold  water  immersion 


rw  (°c) 

20 

28 

Group  Classification 

LliF 

MiiF 

LbF 

MBF 

BMR  (W/m2)* 

49.4 

47.2 

50.5 

42.5 

Tre  (°C)* 

37.48 

37.48 

37.32 

37.52 

^air-neutral 

28.60 

28.85 

28.25 

29.80 

Tbi  (°C) 

37.26 

37.25 

37.10 

37.32 

Tskss(°C)* 

21.0 

21.0 

28.4 

28.4 

1/ 


•HI* 

“7*  I 


Table  2: 

Model  prediction  after  1  h  of  water  immersion 

Tw  (°C) 

20 

23 

Group  Classification 

L13F 

VU3F 

L13F 

MBF 

M  (W/m2) 

221.1 

154.0 

94.3 

64.4 

C  (W/m2) 

214.9 

166.0 

99.7 

32.3 

Hr  (W/m2) 

17.6 

12.9 

6.4 

5.1 

Tre  (°C) 

36.54 

36.61 

36.93 

37.00 

k  (W/m2/°C) 

13.33 

10.63 

11.69 

9  57 

hc  (W/m2/°C) 

214.9 

166.0 

249.3 

205.3 

Tre*  (oc/h) 

-0.37 

-1.13 

-0.36 

-0.64 

V  (oc/h) 

% 

-0.36 

-0.72 

-0.32 

-0.62 

V*(°c/v, 

-0.32 

-0.67 

-0.32 

-0.62 

•slope  of  temperature  against  time 
••calculated  us»n£  Eq.  v. 
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used  In  t to  thermoregulatory  model.  £*ch  body  Moment  la  coo^>oaed  oi 
(our  concentric  anmiiar  compartment*,  the  toed  modelled  a*  a  ipton# 
end  tto  other*  a*  cylinder*.  Length  oi  ths  cylinder*  era  given  to  cm 
U&X  The  central  blood  compartment  U  kxated  «tthto  tto  tnjnfc 
•eg  ment.  Tto  number*  In  perentheeU  repreeant*  t to  outer  radii  (cm)  of 
tto  modal  mg  menu  to r  tto  Uif  group  uned  to  due  (tody. 
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wUe  rats  UA)  plotted  nalnjt  time  (or  the  MBF  group  (n«6) 
in  20°  water.  The  darted  line  rtows  the  prediction  according 
Ttgomerf  controller  lot  rti  rering  (21). 


predicted 


Pkj,  T.  Predicted  meen  body  temperature  (Tj,),  reetei  temperature  (Tre), 
ddjrt  t*rfiptT»tijre  (T^),  metibolic  rate  (J4),  and  convective 
Or and**  (O,  tor  Ove  UftP  group  imrrvwvetl  in  20°C  water. 


APPENDIX  A;  Convective  Heat  Transfer  Coefficient 


The  convective  heat  transfer  coefficient  for  flow  across  spherical  and 
cylindrical  segments  is  determined  by  (21) 

hc  =  Kw  Nu/d  (Al) 

where  Kw  is  the  thermal  conductivity  of  water,  Nu  is  the  Nusselt  number,  and  d 
is  the  segment  diameter.  Convective  heat  transfer  involves  both  forced  and  free 
convection. 

The  Nusselt  number  for  forced  convection  is  determined  by  (7) 

Nuf0  =  0.66  Re1/2  Pr1/3,  (A2) 

where  Re  is  the  Reynolds  number  and  Pr  is  the  Prandtl  number.  The  Reynolds 
number  is  determined  by 

Re  =  Vw  d/  v,  (A3) 

where  Vw  is  the  water  velocity  and  v  is  the  kinematic  viscosity  of  water.  The 
Prandtl  number  is  determined  by 

Pr  =  v/D,  (A4) 

where  D  is  the  molecular  diffusivity  of  water. 

The  Nusselt  number  for  free  convection  is  determined  by  (7) 

Nufr  =  0.54  (PrGr)l/*,  (A5) 

where  Gr  is  the  Grashof  number  determined  by 

Gr  ■■=  BgD3  (Tsk  -  Tw)/v2,  (A6) 

and  where  0  is  the  coefficient  of  thermal  expansion  of  water  and  g  is  the 
acceleration  due  to  gravity. 

If  the  ratio  Gr/Re2  is  small,  then  forced  convection  dominates,  otherwise 
free  convection  dominates  (7).  When  the  ratio  is  near  unity,  it  is  assumed  that 
the  two  terms  are  additive. 
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APPENDIX  A:  AASIC  Statements  Implementing  Thermoregulatory  Model  and  A 
Sample  Simulation 

The  source  program  typed  out  by  the  computer  is  listed  below.  The  first 
section  (lines  10  to  21)  lists  all  the  references  used  to  develop  the  program.  This 
is  followed  by  the  identification  of  body  segments  and  compartments  (lines  29  to 
38)  and  a  glossary  of  terms  not  identified  elsewhere  in  the  program  (lines  30  to 
88). 

Next  begins  the  description  of  the  CONTROLLED  SYSTEM  (lines  100  to 
1260).  Here  given  inputs  of  subject  height,  weight,  percent  body  fat,  and  oasal 
metabolic  rate,  an  extensive  physical  and  physiological  description  of  the  model 
compartments  is  computed.  Compartment  weight,  dimension,  thermal 
capacitance,  basal  heat  production,  basal  blood  flow,  and  thermal  conductance 
are  given  m  order  of  head,  trunk,  arms,  hands,  legs,  and  feet  from  top  to  bottom, 
and  core,  muscle,  fat,  and  skm  from  left  to  right-  An  example  follows  the  source 
listing. 

The  nex;  section  describes  the  CONTROLLING  SYSTEM  (lines  2090  to 
2/20).  Estimates  of  the  otf  temperature  and  central  blood  temperature  tor 
condition  of  thermal  neutrality  are  entered;  if  this  condition  is  not  satisfied  (to 
be  determined  by  the  user  later),  then  the  program  is  re-run  with  new  estimates. 
The  water  temperature,  steady  state  skin  temperature,  net  amount  of  internal 
power  produced  to  exercise,  time  step  for  printout  of  results,  and  the  total 
simulation  lime  are  entered.  Finally,  live*  values  of  the  control  coefficients,  CO 
through  C /»  arc  entered^  CO  specifies  the  half-time  for  the  onset  of  limb 
shivering;  Cl  through  C6  are  the  parameters  that  define  the  shivering  response 
(see  tup  6  m  text);  and  CZ  determines  the  extent  of  counter -current  heat 
exchange.  In  addition  to  calculating  the  neutral  temperature  distribution  from 
which  the  model's  set-point  values  are  assigned  ia  printout  is  optional),  this 


section  specifies  the  skin  thermoreceptor  inputs  and  effector  outputs  and  the 
distribution  factors  of  heat  production  for  muscle  due  to  exercise  and  shivering. 

Simulation  of  the  exposure  begins  on  line  2730.  Warm  and  cold  signals  are 
established  by  comparing  a  compartment's  set-point  value  to  its  current 
temperature.  Initial  signals  upon  immersion  in  cold  water  arise  predominantly 
from  the  skin  compartments  and  to  a  much  lesser  extent  from  core  temperatures 
affected  by  the  respired  heat  loss.  The  efferent  outflow  or  amount  of  shivering 
and  vasoconstriction  arc  determined  in  lines  2930  and  29S0.  The  efferent  output 
or  metabolic  rate,  blood  f  ow,  and  respired  evaporative  and  convective  heat  loss 
are  determined  in  lines  3010  to  3270.  Convective  heat  transfer  coefficients  for 
the  initial  transient  cooling  of  the  skin  are  determined  in  lines  3310  to  35SO. 
Heat  flow  values  for  ail  compartments  including  the  central  blood  arc 
determined  in  lines  36 IQ  to  3790.  The  optional  countercurrent  heat  exchange  is 
declared  by  specifying  C7  as  non-set  s>;  The  integration  step  is  then  determined 
in  lines  3Si0  to  3871  under  the  constraint  that  the  change  in  temperature  of  any 
compartment  cannot  exceed  Q.i^C  during  that  step,  before  the  simulation 
continues,  a  screen  printout  of  the  predicted  initial  trunk  core  temperature  (Tre) 
and  heat  storage  of  the  central  blood  (Fblood)  is  displayed,  if  Tfe  does  not 
match  the  initial  measured  rectal  temperature  or  if  |Fblood| exceed*  0»l®C/min, 
then  enter  new  Estimates  of  f neutral  and  Tblood.  Otherwise,  the  program 
calculates  *  -  new  cotnpaf  tmtm  temperatures  and  the  simulation  continues  until 
the  exposure  over.  A  printout  of  the  model  temperature  distribution  is 
option^;. 

An  example  nan  follows,  this  simulation  was  for  the  W3F  group  immersed 
m  water  (see  Pigs.  2,  $  and  7),  Following  the  description  of  the 

CONTROLLED  SViTLU,  the  neutral  temperature  distribution  is  given  in  the 
compartment*!  order  described  earlier  where  tl,  T2,  O  and  D  refer  to  the 

n 


cores  muscle,  fat  a, id  skin,  respectively.  Following  this  are  ihe  results  for  the 
simulated  immersion  given  every  6  min  up  to  1  h.  TIME  is  the  time  in  min,  TB  is 
the  mean  body  temperature  in  °C,  MR  is  the  metabolic  fate  in  W,  EV  is  the  total 
evaporative  heat  loss  in  W,  TBL  is  the  central  blood  temperature  in  °C,  TCR  is 
the  trunk  core  (rectal)  temperature  in  °C,  THD  is  the  head  core  (hypothalmus) 
temperature  in  °C,  and  HFSK  is  the  convective  heat  loss  in  W/m^  (values  greate- 
than  999  are  indicated  as  999). 
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I  REM  BASIC  listing  of  RWBC  (Thermoregulatory  model  for  Whole  Body  Cooling  of 
nude  subject  immersed  in  cold  water) 

9  REM 

10  REM  References  used  in  the  listina 

II  REM 

IE:  REM  Ref  1  Montgomery  LD.  Annals  Biomed  Eng  v£  1974  o  19-46 

13  REM  Ref  £  Stolwijk  JAJ.  Mathematical  Model  of  Thermorepul  at  ion.  Int 

Physiological  and  Behavioral  Temoerature  Regulation  197®,  o703-72l 

14  REM  Ref  3  Stolwijk  JAJ,  Hardy  JD.  Control  of  Body  Temperature.  In:  Handbook 

of  Physiol gy  1977,  o45— 68 

15  REM  Ref  4  Bullard  RW,  Raoo  GM.  Aerospace  Med  v4i  197®  e  1269-1277 

16  REM  Ref  5  Gagge  AP,  Nishi  Y,  Gonzalez  RR.  Standard  Effective  Temperature  - 

A  Single  Temoerature  Sensation  and  Thermal  Discomfort 

17  REM  Ref  6  Witherspoon  JM,  Goldman  RF,  Breckenridge  JR.  J  de  Physiologic  v63 

1971  0459-462 

16  REM  Ref  7  Stolwijk  JAJ,  Hardy  JD.  Pfluoers  Archiv  v291  1966  pl£9-162 
19  REM  Ref  6  Campbell  GS.  An  Introduction  to  Environmental  Biophysics. 
Sorinqt-Verlag,  1977 

£0  REM  Ref  9  Sekins  KM,  Emery  AF.  Thermal  Science  for  Physical  Medicine.  In  : 

Therapeutic  Heat  and  Cold.  Baltimore.  Williams  and  Wilkins  1982 
£1  REM  Ref  10  Ruch  TC,  Patton  HD.  Physiology  and  Biophysics.  W.  B.  Saunders.  19£'j 
28  REM 

£9  REM  Identification  of  body  segments  and  compartments 
30  REM 

3£  REM  I  refers  to  dodv  segments  as  follows 

34  REM  1  =  head  £  =  trunk  3  =  arms  4  =  hands  5  =  legs  6  =  feet 
36  REM  N  refers  to  segment  comoartrnents  as  follows 
38  REM  1  =  core  £  =  muscle  3  =  fat  4  =  skin 
40  REM 

50  REM  Glossary  of  terms 

51  REM 

5£  REM  FBF  =  fractional  body  fat 

53  REM  SG  =  body  specific  gravity 

54  REM  AT  =  weight  (kg)  of  adipose  tissue 

55  REM  NAT  =  weight  (kg)  of  non-ad i pose  tissue 

56  REM  CB  ~  body  heat  caoacity  <kcal/C> 

57  REM  MR  =  metabolic  rate  (kcal/h) 

58  REM  QSF  (WSR)  =  basal  metabolism  (weight)  of  the  skin  and  fat 

59  REM  QM  (WM)  =  basal  metabolism  (weight)  of  the  muscle 

60  REM  QC  <WC)  -  basal  metabolism  (weight)  of  the  skeleton  &  connective  tissue 

61  REM  TR  =  thermal  resistance  <C*h/kcal) 

62  REM  K  tissue  thermal  conductivity  (kcal/h/m/C) 

63  REM  Tneutral  *  air  temperature  repuired  for  thermal  neutrality 

64  REM  Tblood  “  central  blood  temoerature  reauired  to  obtain  correct  Tre 

65  REM  Note  that  Tneutral  and  Tblood  are  estimated  by  trial  and  error  until  hea 

storage  of  blood  =  0  and  7(11)  =  initial  measured  rectal  temoerature 
£6  REM  SSTSK  =  steady  state  skin  temperature 

67  REM  WORK I  «  net  amount  of  internal  power  produced  by  exercise  < kcal/h) 

60  REM  DP'RT  =  time  increment  (min)  for  printout  of  prediction 

69  REM  TMAX  “  total  simulation  time  (min) 

70  REM  TIM  =  time  (h>;  TIME  ■  time  (min) 

7 1  REM  7PRT  **  time  of  or  int  out  <h) 

72  REM  PAIR  =■  vaoour  pressure  (rnrnHg) 

73  REM  TB  =  mean  weighted  body  temoerature  (C) 
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74  REM  EV  ■=  evaoorative  heat  lose  <kcal/h) 

75  REM  TBl  *»  central  blood  temoerature  <C> 

76  REM  TCR  »  trunk  core  temoerature  <C) 

77  REM  THD  =  head  core  temoerature  <C> 

70  REM  HFBK  «  convective  neat  loss  to  water  ( kcal /h/m^S) 

79  REM  DIFF  =  difference  between  compartment  temoerature  and  its  set-ooint 
00  REM  WARM  (COLD)  «  warm  (cold)  signal 
01  REM  Q  ==  metabolic  rate  (kcal/h) 

82  REM  BF  *s  blood  flow  (1/h) 

03  REM  E  =  svaoorative  heat  loss  <kcal/h) 

84  REM  BC  --  conductive  heat  exchange  between  comoartment  and  central  blood 

05  REM  TD  «  conductive  heat  exchange  between  compartments 

BE/  REM  HF  =  heat  storage  of  comoartment 

87  REM  F  =  rate  of  change  of  compartment  temoerature 

38  REM  X  ~  factor  to  introduce  limb  shiverinn  exoonent ial ly  with  half-time  CC 

90  REM 

91  LPRINT  "RESTING  WHOLE  BODY  COOLING"  :  LPRINT  "  "  :  LPR1NT  "  " 

92  REM 

100  REM  Description  of  the  CONTROLLED  SYSTEM 
110  REM 

120  REM  Calculate  Body  Weight  Distribution  from  Pierson  and  Eagle  (1915), 
see  Ref  1  o24 

130  REM 

140  INPUT  "Subject  height  (cm),  weight  (kg),  and  body  fat  OO";  HT, WT, PBF 

150  FBF  =  PBF/ 100 

160  SG  --=  5.  548001/ (FBF  +  5.044) 

170  IF  FBF  <  0  THEN  PRINT  "ERROR  LINE  140":  END 
180  AT  =  FBF*WT 
190  NAT  =  WT  -  AT 
200  REM 

210  REM  Calculate  Surface  Area  <SA)  from  DuBois  ana  DuBois  (1915),  see  Ref  ! 
and  Assign  Segment  Lengths  <L>,  see  Ref  £  d708 

220  REM 
225  DIM  L (7) 

230  SA  =  .  0071 84*WT  ''.  425*HT-''.  725 

240  L (2)  =  .6 

250  L (3)  -  1. 12 

260  L (4)  ■  .96 

270  L  <  5)  =  1.6 

280  _< 6)  =1.25 

290  REM 

300  REM  Table  of  Relative  Weight  Distribution  for  Model  Compartments  (CSWT) , 


see  Ref 

l 

025 

310 

REM 

320 

DIM  CSWT  <  7 , 

lit 

3) 

340 

cswr  (1,  l,  1 ) 

= 

.  0262:32 

350 

CSWT  ( 1 ,  1,2) 

* 

. 02518 

360 

CSWT  (1,2,  1  ) 

— 

. 00588 

370 

CSWT (2, 1,1) 

» 

. 18704 

380 

CSWT  <2,  1,2) 

at 

. 038832 

390 

CSWT (2, 2, 1) 

a 

.  2834 

400 

CSWT (3, 1,1) 

■* 

. 01)764 

410 

CSWT (3, 1, 2) 

= 

. 023756 

420 

CSWT <3, 2, i) 

— 

. 05328 
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<4  30 

CSWT  <4,  1,1) 

xr-  .  0004704 

440 

CSWT  <4,  1, 2) 

=  . 003646 

4  50 

CSWT <4, 2, 1) 

=  .001188 

460 

CSWT  <5, 1.1) 

=  . 030352 

470 

CSWT  <S,  1,2) 

=  .079168 

480 

CSWT (5, 2, 1 ) 

=  .  1 6 1 

490 

CSWT  <6,  1,1) 

=  .0009412 

500 

CSWT  <6,  1,2) 

=■  .  005664 

510 

CSWT <6, 2, 1 ) 

=  . 00 1 1 83 

530 

CSWM1,  3,  1) 

=•  .  0333 

540 

CSWT <1, A, 1) 

-  .  004 £3 

550 

CSWT  <2, 3,  1 ) 

-  . 6333 

560 

CSWT  <2, 4,  1) 

=  .0213 

570 

CSWT  (3,3,  1 ) 

■-  .  0867 

580 

CSW  r  <3, 4,  1) 

■=  .00764 

590 

CSWT  <4,  3,  1 ) 

-  .01333 

600 

CSWT  <4, 4,  1 ) 

=  .00294 

610 

CSWT (5, 3, 1 > 

—  •  d  jw  Sj  v  j 

620 

CSWT  <5, 4,  1 ) 

=  .01894 

630 

CSWT  <6,  3,  1) 

-  .  08 

640 

CSWT <6, 4, 1) 

®  .00376 

650  REM 

660  REM  Calculate  Compartment  Weights,  SWT(kg>,  and  Thermal  Capacitance 
C(kcal/C),  see  Ref  1  p26 

670  REM 

660  DIM  SWT(70),C<70) 

685  CB  =  0 

690  FOR  I  *  1  TO  6 

700  FOR  N  =  1  TO  4 

710  SUIT  <  1 0* I +N—  1 0 )  (CSWT  <  I ,  M,  1  )  +  CSWT  ( I ,  N,  2)  )  *NflT 

720  C  ( 1 0*  I +  N- 10)  =  (.  9*CSWT  <  I,  N,  1  >  +  .  5*CSWT  ( I ,  *v,  2)  )  *N«T 

730  NEXT  N 

740  SWT(10*I-7)  =  CSWT (1,3,  i >  *AT 

750  C<  10*1-7)  =  .  £*SW?  ( 1 0* I  —7 ) 

755  CB  «  CB  +  C  < 10*1—9)  *  C(l0*l-8>  +  CU0*I-7)  +•  C(10*I-&> 

760  NEXT  I 

770  REM  to  account  for  the  thermal  caoacitance  of  blood  in  trunk  core, 
see  Ref  2  o708 
780  C<11)  *»  C(ll)  -  2.  25 
790  C  <61  >  -  2.25 

800  REM  Calculate  thermal  capacitance  of  immersed  skin  <CS) 

810  CS  =  0 

620  FOR  I  =  1  TO  6  :  CS  =  CS  +  C< 10*1-6)  s  NEXT  I 
825  REM 

S30  REM  Calculate  Basal  Heat  Production,  QEMkcal/h),  see  Ref  1  p27 
640  REM 

850  INPUT  "Resting  metabolic  rate  (kcal /h/nr'l£)  " ;  BMR 
860  MR  =  BMR*SR 

670  QSF  *  . 3* (. 05SB2*NfiT  +  QT) 
fl60  QM  =  .  1S*MR  -  QSF 
690  DC  ~  .  1 *MR 

900  WC  ~  0  :  WM  =:  0  :  WSF  *=  0 

910  FOR  I  -  1  TO  6 

930  WC  «  WC  +  SWT < 10*1-9) 


Va 1 ues , 
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940  WM  =  WM  +  SWT < 10* I -8) 

960  WSF  =  WSc  +  SWT ( 10*1—7)  +  SWT (10*1-6) 

970  NEXT  I 
980  DIM  QB  < 70 ) 

1000  FOR  1=1  TO  8 

1010  QB ( 10*1—9)  =  SWT < 10*1-9) *QC/WC 
1020  QB  (10*1-8)  =  SWT < 10*1-8) *QM/WM 
1040  QB ( 10*1—7)  =  SWT < 10*1-7) *QSF/WSF 
1050  QB ( 10*2— &)  =  SWT ( 10*1 -8) *QSF/WSF 
1080  NEXT  I 

.1070  REM  to  account  for  extra  values  for  head  ana  trunk  cores,  see  Ref  1  p£7 
1090  QB  (  1 )  =  DBC)  +  . 16*MR 
1 1 00  QB  (.11)  =  QB  <  U  >  *  .  S6*MR 
1115  QC  =  . 82*MR 
1120  REM 

1130  REM  Calculate  Basal  Blood  Flow,  BFB<l/h),  see  Ref  1  p£8 

1140  REM 

1150  DIM  BFB (70) 

1160  FOR  1=1  TO  6  :  FOR  N  =  1  TO  4 
1170  BFB ( 10*I+N-10)  =  1. £*QB( 10*I+N-10) 

1180  NEXT  iM  :  NEXT  I 
1 £00  BFB ( 1 )  =  45 
1210  BFB (11)  =  £10 
1230  BFB (4)  =  5. 34*SWT  < 4 ) 

1240  BFB (14)  =  1 „ 56*BWT (14) 

1250  BFB (£4)  =  1.04*SWT(24) 

1260  BFB (34)  =  10. 5*SWT<34) 

1270  BCB<44>  =  2.  38*SWT  (44.' 

1280  BFB (54)  =  12. 5*SWT (54) 

1 300  REM 

1310  REM  Calculate  Compart merrc a  1  Volume,  V(rn"-3),  Radii,  R(rn>,  Interfacial  Areas, 
fi<nv''£),  and  Thermal  Conductances,  !'C (kca  1 /C/h ) ,  see  Ref  £  d707  and 
Ref  9  Q84&94 

1315  REM  ftssurnes  volume  (riM'S)  **  weight  <kg)/1000 
1320  REM 

1330  DIM  V  (7)  ,  R  (70)  ,  A(70)  ,  RCM  (70)  ,  TC  (70)  ,  TR  (70.)  ,  K  (70) 

1332  FOR  I  =  1  TO  6 

1334  K( 10*1-9)  =  .  36  :  K< 10*1-8)  =  .2394 
1336  K  < 10*1—7)  =  .1634  :  K ( 10*1-6) '  «  .288 
1338  NEXT  1 

1340  K ( 1 2 )  =  .2988  :  K<£2>  =  .2988  :  K(4£)  =  .2938 
1345  PI  =  3. 14159 
1350  FOR  I  =  1  TO  6 
1360  FOR  N  =  1  TO  4 

1370  V<I>  =  V(I)  +  SWT  < 1 0*1 +N— 10) / 1 000 
1380  NEXT  N 
1390  NEXT  I 

1400  REM  head  segment  as  sohere 
1410  VOL  =  V ( i ) 

1420  R  (  4  )  =  (.  75*V0L/PI  )•''<! /3) 

1425  RCM  <  4)  =  <R<4)~3  -  . 000375*SWT (4 ) /PI ) ~ ( 1 /3> 

1430  ft  <  4 )  »  4*PI *R  < 4 )  ',%£ 

1440  FOR  N  ~  1  TO  3 

1450  VOc  =  VOL  -  SWT (5-N) / 1000 
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A  460  R(4-N)  =  (.  75*V0L/PI )  (  1/3) 

1465  RCM (4— N)  «*  (R(4-N)\3  -  . 080375*SWT  <4-N> /PI  >  ^  U/3) 

147©  ft  <  4— N)  =  4*PI*R<4-N)"'£ 

1480  NEXT  N 

15£0  FOR  N  =  1  TO  3 

1530  T  R  <  N )  =  <1/RCM(N>  -  1/R (N> ) / (4*PI*K (N) >  +  (1/R(N)  -  1  /  RCM  <N+ 1 )  )/ ( 4*PI*k  ( N+ 
i  >  ) 

1540  TC(N)  =  1/TR(N) 

1550  NEXT  N 

1560  REN  remaining  seaments  as  cylinders 
1570  FOR  I  «  £  TO’ 6 
1580  VOL  =  V(I) 

1590  R( 10*1-6)  =  (VOL/ <PJ*L ( I) ) ) A. 5 

1595  RCM  ( 10*1-6)  *»  (R  (  10*1-6)  ■''S  -  .  0005*SWT  ( 10*  I -6)  /  < PI  *L  ( I  >  >  > '•.  5 
1600  ft ( 10*1-6)  »  £*PI*R ( 10*1-6) *L (I ) 

1610  FOR  N  =  1  TO  3 

16S0  VOL  »  VOL  -  SWT < 10*1 -N-5) / 1000 
1630  R(10*I-N-6)  »  <VOL/(PI*L(I) ) )A.S 

1635  RCM  <10*1  — N— 6)  =  (  R  ( 10*1— N-6)  •*'£  -  .  0005*SWT  ( 10*I-N-6)  /  <PI*L  ( I )  )  ) 's.  5 

1640  ft ( 10*1— N— 6)  =  £*PI*R ( 10*1 -N-6) *L ( I ) 

1650  NEXT  N 

1700  FOR  N  =  1  TO  3 

1705  J  =  10*I+N— 10 

1710  TR(J)  =  LOG(R(J) /RCM(J) )/ (S*PI*K ( J) *L < I ) )  +  LOS ( RCM < J+l ) /R ( J) ) / (£*PI*K < J+l 
)  *L  <  I )  ) 


17E0 

TC  ( J  > 

=  1/TR(J) 

1730 

NEXT  N 

1740 

NEXT  I 

1750 

REM 

1760 

REM  End 

of  Descriotion 

1770 

REM 

1780 

INPUT  " 

Enter  1  for  model  descriotion,  else  0";  CODE 

1790 

IF  CODE 

=  0  THEN 

GOTO  £100 

1795 

LPRINT 

II  ll 

1800 

LPRINT 

"H7  ="  HT 

"  WT  ="  WT  " 

SG  ="  SG  "  FBF  «"  FBF  "  Sft  =»"  Sft 

1810 

LPRINT 

"QC  ="  QC 

"  QM  QM  " 

QSF  ="  QSF 

1615 

LPRINT 

"WC  ="  WC 

••  WM  ="  WM  " 

WSF  ="  WSF 

1 8£0 

LPRINT 

t»  U 

1S£5 

LPRINT 

"  Core 

Muscle  Fat 

Skin" 

1830 

LPRINT 

"  "  :  LPRINT  "Weight  SWT 

(kg)"  :  LPRINT  »  » 

1840 

FOR  I  = 

1  TO  6 

1850  LPRINT  USING  "##.###  "•  SWT ( 10*1-9)  ,  SWT < 10*1 -8) , SWT ( 10*1-7) , SWT (10* 1—6) 

I860  NEXT  I 

1870  LPRINT  M  "  :  LPRINT  "Radius  R  (cm)"  t  LPRINT  "  " 

1880  FOR  I  =  1  TO  6 

1890  LPRINT  USING  »##.  ##  ";  100*R  (10*1-9)  ,  1 00*  R  (10*1-8)  ,  100*R  ( 10*1-7)  ,  .100*  R  (  1 

0*Ir6) 

1900  NEXT  1 

1910  LPRINT  "  »  :  LPRINT  "Thermal  Caoacitance  C  (kcal/C)"  :  LPRINT  "  " 

1920  FOR  I  =  1  TO  6 

1930  LPRINT  USING  »*)#.###  ";  C  <  10*1  -9) ,  C  ( 1 0*1-8) ,  C  <  1 0*1-7)  ,  C  ( 10*1-6) 

1940  NEXT  I 

1950  LPRINT  "  "  j  LPRINT  "Basal  Heat  Production  OB  (kcal/n) "  :  LPRINT  "  " 

I960  FOR  I  -  1  TO  6 
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1970  LPRINT  USING  "##.###  QB  ( 10*1  -9) ,  QB  <  10*  I -8) ,  QB  <  1 0#  1-7 ) ,  OB  (  i  0*1-6) 

1 980  NEXT  I 

1990  LPRINT  "  "  s  LPRINT  "Basal  Blood  Flow  BFB  <l/h)"  s  LPRINT  ”  " 

£000  FOR  1  ■  i  TO  6 

£010  LPRINT  USING  "###.##  BFB < 10*1-9)  ,  BFB ( 10*1-8)  ,  BFB < 10*1-7)  ,  BFB ( 10*1-6) 

£020  NEXT  I 

£030  LPRINT  "  •’  s  LPRINT  "Thermal  Conductance  TC  <kcal/C/h)"  :  LPRINT  "  " 

£040  FOR  I  *  l  TO  6 

£050  LPRINT  USING  •'##.##  "  ;  TC  <  1 0*1-9)  ,  TC  ( 1  0*1-8)  ,  TC  ( 10*1—7)  ,  TC  ( 10*1—6) 

£060  NEXT  I  • <  • 

£065  OB  =  0(4)  +  0(14)  +  A(£4)  +  0(34)  +  A<44)  +  0(54) 

2070  LPRINT  "  "  :  LPRINT  "Segment  Surface  Area  0  (m*m) "  :  LPRINT  "  " 

£080  LPRINT  "01  ="  0(4)  "  ft£  ="  0(14)  "  03  =  "  0(£4)  "  04  ="  0(34) 

£081  LPRINT  "05  *"  0(44)  "  06  »"  0(54)  "  OB  ="  OB  :  LPRINT  "  " 

£090  REN 

£100  REN  Description  of  the  CONTROLLING  SYSTEM 
£105  REM 

£110  REN  Initial  conditions,  assumes  relative  humidity  (RH)  *=  1,  and  air 
temperature  (TOIR)  =  water  temperature  (TWOT) 

£ 1 £0  REM 

£125  BIN  EB  (70)  ,  HSS<7)  ,  T  (70)  ,  TSET  (70) ,  TF(il,  1 1 )  ,  EF  ( 1 1 )  ,  W0RKM(7)  ,  CHILM(7> 

£1£6  DIN  SKI NR (7) , SKINS (7) , SKINV (7)  ,  SKINC (7) 

£1£7  DIN  WORN (70) , COLD  <70) , DIFF (70) 

£l£8  DIM  Q (70)  ,  BF (70) , E(70) , TBL (70) 

£1£9  DIN  BC<70) , TD<70) , HF (70) , F (70) 

£130  INPUT  "Tneutral  Tblood  Twater  ssTSK  WORKI  DPR I NT  TMOX";  TO,TB,TWAT,S 
TSK, WORKI, DPRT,  TMOX 

£135  DPRT  =  BPRT/60  :  TMOX  *  TMOX/60 
£140  TOIR  =  TWOT 

£150  i-'O  I R  =  EXP(18.  6686  -  4030.  183/ (TOI R  +  £35)) 

£160  TIM  =  0  :  TPRT  =  0 

£165  LPRINT  "  "  :  LPRINT  "Tneutral-air  ="TP"C  Tblooc-  ="TB"C  Twater  «"TW0T"C 
:  LPRINT  "  » 

£170  REM  Input  control  constants;  C0  see  £905;;  Cl  to  C6  see  £930;  C7  see  3637 
£180  INPUT  "Enter  control  constants  C0  to  C7" ;  C0, C 1 . C£, C3, C4. C5. C6, 07 
£185  REM 

£190  LPRINT  "  "  :  LPRINT  "C0  ="C0"  Cl  ="C1"  C£  ~"C£"  C3  ~"C3"  C4  =  »C4"  CL 
"C5"  C6  ="C6"  C7  ="C7  :  LPRINT  "  " 

££00  REM 

££30  REM  Calculate  initial  temperature  distribution,  T(N),  assuming  subject  i 
thermally  neutral  for  given  Tneutral-air  (RH  =  .5)  and  Tblood 

££3£  REM 

££34  REM  Table  of  basal  evaporative  rates,  EB(kcal/h),  see  lines  3160-3£30  ana 
Ref  £  p  708  and  Ref  5  o  £46 
££38  PA  =  . 5*EXP( 18. 6686  -  4030. 183/ (TA  +  £35)) 

££40  ERES  =  . 00197B*MR*<44  -  PA) 

££4£  CRES  *  . 00103£*MR* (34  -  TO) 

££44  EB (4)  =  .6120001  :  EB(14)  =  3.27  :  EB(£4)  =1.165 
£245  EB ( 34 )  =  . 43£  :  EB<44)  =  £.  98  :  EB<54)  =  .6 
£246  EB  < 1 >  =  4.5  :  EB ( 1 1 )  =  ERES  +  CRES 
££50  REM 

££5£  REN  Table  of  heat  transfer  coefficients  from  BKin  to  air,  see  Ref  3  p  5r* 
££54  REN 

£256  HSS(l)  =  6.7140(4)  s  HSS(£)  »  5.93*0(14) 
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££56  HSS  (3)  =  5. 4£#A<24>  s  hSS<4)  *  6. 53*A(34> 

££60  HSS (5)  =  5. 4£*A < 44 )  :  *55(6)  =  6. 88*fi(54> 

££•90  REM 

££9£  REM  Calculate  comoart mental  temoerature  and  matrix  coefficients  for  the 
4x4  reoresentation  of  the  non-homogeneous  system  of  eouat  ions  for 
determining  the  neutral  temperature  distribution 

££93  REM 

£300  FOR  I  =  1  TO  6 

£302  TF  (1,1)  =  . 9*BFB( 10*1-9)  +  ~C  (  .10*1-9) 

£304  TF (1,2)  =  -TC (10*1-9) 

£306  EF ( 1 )  =  QB (10*1-9)  -  EB < 10*1-9)  +  . 9*BFB < 10*1-9) *TB 
£308  TF (4,3)  =  -TC(10*I-7> 

£310  TF  (4,  4 )  *  .  9*BFBU0*I-6)  +  TC(10*I-7)  +  HSS  ( I ) 

£3 l £  EF ( 4 )  =  QB ( 10*1—6)  -  EB(l0*l-6>  +  . 9*BFB ( 10*1-6) *TB  +  HSS < 1 > *TA 

£314  FOR  N  =  £  TO  3 

£316  TF  (I'm,  N— 1 )  =  -TC  ( 10*I+N-1 1 ) 

£316  TF<N, N)  =  . 9*BFB < 10*I+N-10)  +  TC < 1 0* I +N-1 1 )  +  TC  < 10*1 +N— 10) 

£3£0  TF (N, N+l )  *  -TC ( 1©*I+N-10) 

£3££  EF  <  N)  ••=  QB  <  10*1  +N-10)  -  EB < 10*I+N~10>  +  . 9*BFB ( 10*I+N-10) *TB 
£3£4  NEXT  N 
£3 30  REM 


£33£  REM  Begin  Gauss  Elimination  method  by  initializing  first  terms  to  unity 
normalizing  accordingly 


c!u)owi 

REM 

£334 

EF  ( 1 )  = 

EF(1)/T 

F  < 

1,  1) 

£336 

TF  < 1 , S) 

=  TF  <  1 , 

£ ) 

/TF ( 1 , 1) 

£336 

EF<4)  = 

EF (4) /T 

=  < 

4,  3) 

£340 

TF  <4, 4) 

=  TF (4, 

4) 

/TP (4, 3) 

£34£ 

FOR  N  = 

£  -rO  3 

£344 

EF(N>  = 

EF (N) /TF 

(N,  N-l ) 

£346 

TF (N,  N+ 

1)  »  TF 

(N 

,  N+l  )  /TF  CM,  N- 

£346 

TF (N, N) 

=  TF  (N 

) /TF (N, N— 1 ) 

£350 

NEXT  N 

£355 

REM 

£360  REM  Begin  elimination  procedure  and  normalize 
£361  REM 


£36£ 

FOR  N  = 

£  TO  3 

£364 

TF  <N,  N) 

-  TF  <N,  N)  -  TF  < N— 1 , N> 

£366 

EF (N)  = 

EF (N)  -  EF(N-i) 

£368 

EF (N)  = 

EF(N)/TF(N,N) 

£370 

TP (N, N+l)  =  TF (N, N+l > /TF <N, N> 

£37£ 

NEXT  N 

£374 

TF (4, 4) 

=  TF  (4,  4)  -  TF  (3,  4) 

£37.6 

EF ( 4)  = 

(EF (4)  -  £F(3> ) /TF (4, 4) 

2376 

REM 

£360 

REM  Begin 

suast it ut ions 

£38.1 

REM 

£38£ 

T< 10*2-6)  =  EF  < 4 ) 

£384 

FOR  N  = 

1  TO  3 

£366  T ( 10*1— N— 6)  =  EF(4-N>  -  TF (4— N,  5~N> #T ( 10*1-M— 5) 
£366  NEXT  N 
£390  NEXT  I 
£400  REM 

£40£  REM  fission  set-point  iemoeratures,  TSE’f  <C) 

40  . 


and 


£403  REM 

5404  T(6l)  =.  TB 

5406  -OR  N  »  1  TO  61  s  TSET(M)  ■  T  < N)  s'  NEXT  N 

£410  INPUT  "Enter  l  for  neutral  temoerature  distribution*  else  0";  CODE 
£41£  IF  CODE  *  0  THEN  GOTO  £530 

£414  LPRINT  "  "  :  LPRINT  "  T1  T£  T3  T4"  s  LPRINT  "  " 

£416  FOR  I  *  1  TO  6 

£416  LPRINT  USING  "##.##  ";  T ( 10*1 -9) , T < 1 0* I -8) , T (1 0*1 -7) , T ( 1 0*1 -6) 

£430  NEXT  I 
£530  REM 

£530  REM  Table  of  skin  thermoreceotor  inouts  and  effector  outputs, 
see  Ref  1  p3£ 

£540  REM 


£560 

3KINRU  ) 

=5 

.0695  :  SKINR  <£) 

«  .4935 

i  SKINR (3)  «  .0686 

£570 

SKINR (4) 

= 

.1845  :  SK I NR  <5 > 

a  . 1505 

:  SKINR  <6> ,  =  .  0334 

£580 

SK 1 NS  < 1 ) 

= 

.081  :  SKINS <£> 

=  .481 

SKINS (3) 

=  .  154 

£590 

SKINS  <4) 

St 

.031  :  SKINS (5) 

=  .  £18 

SKINS (6) 

a  .  035 

£600 

SKINV<1) 

= 

. 13£  t  SKINV <£) 

a  .  3££ 

SKINV  (3) 

«  9. 50000 1E“0£ 

£610 

SK INV  < 4) 

= 

.121  :  SKINV (5) 

=  .£3 

SKINV (6) 

a  .  1 

£6£0 

SKINC  < 1 ) 

= 

.05  :  SKINC <£)  = 

.  15  : 

SKINC  <3)  = 

.05 

£630 

SK I NC ( 4 ) 

= 

.35  :  SKINC (5)  ■ 

.  05  : 

SKINC  <6>  =■ 

.  35 

£640  REM 

£650  REM  Table  of  distribution  factors  of  heat  oroauction  for  muscle  due  to 
exercise  (WORKM)  and  to  shiverinp  (CHILM),  see  Ref  1  a33 

£660  REM 

£680  WORKM <1)  =  0  s  WORKM <£)  =  .3  :  WORKM (3)  •=  .08 

£690  WORKM (4)  »  .01  s  WORKM (5)  ®  .6  :  WORKM (6)  =  .01 

£700  CHILM(l)  -  .0Z  :  CHILM <£l  *  .8499999  :  CHILM<3>  =  .05 

£710  CHIi_M<4>  =  0  :  CHILM  (5)  =  .  07  :  CHILM  <6>  =  0 

£715  REM 

£7£0  INPUT  "Enter  1  for  temoerature  a i st r i but  ion,  else  0";  CODE 
£7£5  REM 

£730  REM  Start  of  simulation 
£740  REM 

£745  LPRINT  "  »  :  LPRINT  "  TIME  TB  MR  EV  TBL  TCR 

TSK  HFSK"  :  LPRINT  "  " 

£750  REM  Establish  thermoreceotor  outout  and  interprate  peripheral  efferents 
see  Ref  1  o  31-3£ 

£760  REM 

£775  INIT  =  1 

£780  FOR  n  =  1  TO  60 

£790  WARM(N)  »  0  :  COLD<N)  =0 

£800  D1FF(N)  =  T(N)  -  TSET (N) 

£810  IF  DIFF(N)  >  0  THEN  WARM(N)  =  BIFF  <  N>  ELSE  COlD(N)  *  -DIFF(N) 

28£0  NEXT  N 

£830  WARMS  =  0  :  COLDS  =  0 
£840  FOR  I  «  1  TO  6 

£850  WARMS  =  WARMS  +  WARM < 10*1-6) *SKINR ( I ) 

£860  COLDS  =  COLDS  •*  COLD  <  10*1-6)  *SKINR  ( I  > 

£870  NEXT  I 
£880  REM 

2890  REM  Determine  efferent  outflow,  see  Ref  1  o  33 
£900  REM 

£905  X  *  EXP  <-. 693*60*TI M/C0) 


'rHD 
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2907  CHILM  (2)  «  .  0499999  -*■  .  12*X  3  CHILM(3>  *>  .  05*  <  i  -  X)  :  CHILMCS)  »  . 07* (1  - 
X) 

£910  SWEPT  =  32*DIFF  < 1 )  +  £9*  (WARMS  -  COLDS}) 

£9£0  DILAT  =*  117*DIFF<1>  +  7. 5*  (WARMS  -  COLDS) 

2930  CHILL  =  SA/  I.  163*  (Ci*COLDS-/'C3/PBF/'C4  +  C2*C0LD  ( 1 )  ■»C0LDS^CS/PBF-''C6> 

2940  STRIC  =  -5*DIFF  < 1 )  +  5* (COLDS  -  WARMS) 

2950  IF  SWEAT  <  0  THEN  SWEAT  =  0 

2960  IF  DILAT  (  0  THEN  DILAT  «  0 

£970  IF  CHILL  <  0  THEN  CHIU-  =  0 

2900  IF  STRIC  <  0  THEM  STttJC  =  0 

3000  REM 

3010  REM  Assign  efferent  out  out,  see  Ref  1  D  34-35 
3020  REM 

3040  FOR  I  =  1  TO  6 

3060  Q< 10*1-9)  =  QB  C 10*1-9) 

3070  0(10*1-6)  =  QB( 10*1-8)  +  WGRKM < I ) *W0RKI  +  CHILM < I ) *CHILL 

3080  BF< 10*1-9)  «  BFB< 10*1-9) 

3090  BF  ( 1 0*  1  —8)  =  BFB  ( 1 0*  I  -8 )  ->  0(10*1-8)  -  DB(10*I-8> 

3110  0(10*1-7)  =  QB ( 10*1—7) 

3120  0(10*1-6)  =  QB ( 10*1-6) 

3130  BF ( 10*1—7)  =  BFB ( 10*1-7) 

3140  BF  ( 10*1-6)  =  (  (BFB  <  10*1-6)  SKINV  ( I )  *DILAT)  /  ( 1  *  SKINC  ( I )  *STRIC)  )  *£•'•  (DIF? 
( 10*1-6) /6> 

3150  NEXT  I 

3160  REM  evaoorative  heat  loss  from  the  head,  see  Ref  2  o  708  &  712 
3165  E(l)  =  EBd) 

3180  REM  res oi red  evaporative  and  convective  heat  loss,  see  Ref  5  p  £46 
3190  ERES  =  . 001978*MR* (44  -  PAIR) 

3200  CRES  =  .  00.1032*MR*(34  -  TAiR) 

3220  E ( 1 1 )  =  ERES  +  CRES 

3230  REM  Calculate  total  metabolic  rate 
3240  MR  a  0 

3250  FOR  I  =  1  TO  6  :  FOR  N  -•=  1  TO  4 
3260  'MR  =  IMR  +  Q(10*I+N-10> 

3270  NEXT  N  :  NEXT  I 
3300  REM 

3310  REM  Determine  heat  transfer  coefficients  from  skin  to  water  HSS  (kcal/h/C) 
for  initial  non— steady  state  coolinp,  see  Ref  6  and  Ref  8  o  65—70 

3320  REM 

3370  FOR  I  =  1  TO  6 

3380  DELT  =  T( 10*1-6)  -  TWAT 

3400  REM  Initial  heat  flow  is  approximated  by  assuming  minimal  forced  conveetior 
-  after  steady  state  is  attained,  heat  flow  is  determined  by  assuming 
zero  heatstorace  for  the  immersed  skin 
3410  VEL  =  18 

3420  REM  linear  aoor  of  kinematic  viscosity  <cm*cm/s>  from  CRC 
3430  VU  -  .0148  -  . 0082£4#TWAT 

3440  REM  Reynolds  No.  where  factor  .36  converts  VU  to  m*m/h 
3450  RE  »  2*R < 10*1-6) *VEL/ (. 3G*VU> 

3460  REM  linear  aoor  of  Pranotl  No.  based  on  9.5  at  10C  and  7.3  at  22C 
3470  PR  =  11.3  -  . 18*TWAT 

3480  REM  Nusselt  No.  for  forced  convection,  see  Ref  8  o  66 
3490  NUF0  =  . 66*REA.  5*PR^(l/3) 

3500  REM  linear  aoar  of  thermal  expansion  of  water  (I/O  from  CRC 
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3310  TE  »  9.  999999E~06*TWAT 

35£0  REM  Grashof  No.  where  factor  10''6  to  convert  R(rn)  to  R(crn)  is  imbedded  in 
coefficient  of  NUFR,  see  Ref  Q  d  65 
3530  GR  «  980*TE*(£*R<10*I-6)  >  ''3*(T<10*I-&>  -  TWAT  >  / VU'-2 
3540  REM  Nussel t  No.  for  free  convection,  see  Ref  0  o  69 
3550  NUFR  =>  17.  1*  (GR*PR)  £5 

3560  REM  heat  transfer  coefficient  assuming  Kw  «  0. 5£  kcal/m/h/C  from  CRC 
3570  HSS(I)  =  . 52*  (NUFO  +  NUFR) *A < 10*1-6) / (£*R ( 10*1-6) ) 

3580  NEXT  I 
3600  REM 

3610  REM  Calculate  heat  flows  HF  (kcal/h),  see  Ref  1  o  36 
36£0  REM 

36£5  REM  Calculate  blood  temperatures  for  optional  counter-current  heat  exchant 
(set  C7  =  0  for  no  effect) 

3630  FOR  1  ®  1  TO  6  :  FOR  N  =  1  TO  4 
3632  TBL(10*I+N-10)  =  T<61) 

3634  NEXT  N  s  NEXT  I 

3635  FOR  1  =  3  TO  6  :  FOR  N  =  1  TO  3  :  J  *  10*I+N-10 
3637  TBL  <  J  >  =  T<J>  +  <T<61)  -  T ( J> > *EXP (-C7*C0UDS) 

3639  NEXT  N  :  NEXT  1 

3650  FOR  I  =  1  TO  6 

3655  FOR  N  =  1  TO  3  :  J  -  10*1+N-10 
3660  BCCJ)  =  . 9*BF(J)*(T(J>  -  TBL(J)> 

3670  TD  < J)  =  TC ( J ) * ( T  <  J)  -  T(J+1)> 

3680  NEXT  N 

3690  BC ( 10*1—6)  =  . 9*BF < 10*1-6) * (T (10*1-6)  -  TBL(10*I-6)) 

3700  TD( 10*1-6)  =  HSS(i)*(T (10*1-6)  -  TWAT) 

3710  HF ( 10*1—9)  =  Q  < 10*1—9)  -  6(10*1-9)  -  BC(10*I-9)  -  TD(10*l-9) 

3720  FOR  N  =  £  TO  4 

3730  HF ( 10*I+N— 10)  =  Q(10*I+N-10)  -  E < 1 0* 1 +n- l 0 )  -  BC < 1 0* 3 +N- 1 0)  +  TD(10*I+N- 
1)  -  TD < 10*I+N— 10) 

3740  NEXT  N 
3750  NEXT  1 
3760  HF (61 )  *  0 

3770  FOR  I  =  1  TO  6  :  FOR  N  =  1  TO  4 
3780  HF (61 )  =  HF (61 )  +  BC < 10*I+N-10) 

3790  NEXT  N  :  NEXT  I 
3800  REM 

3810  REM  Determine  ootimum  integration  steo,  BT  (change  in  T  cannot  >  0.  1  C) 
3820  REM. 

3830  DT  =  TPRT  -  TIM 

3840  FOR  I  =  1  TO  6  :  FOR  N  =  1  TO  3 

3850  F<10*I+N— 10)  =  HF ( 10*I+N-10) /C < 10*I+N-10) 

3860  IF  .  1/ABS (F ( 10*I+N-10) )  (  DT  THEN  DT  =  .  1 /OBS (F  (  j  0*I+N— 10) ) 

3862  NEXT  N 

3864  F  < 10*1—6)  =  HF  < 10*1—6) /C ( 10*1—6) 

3866  IF  (T( 10*1-6)  -  TWAT)  (  2  THEN  GOTO  3870 

3868  IF  . l/ABS(F(10*I-6) )  (  DT  THEN  DT  -  . 1 /ABS <F < 10*1-6) > 

3870  NEXT  I 

3871  F (61 )  =  HF  <61)/C(6l) 

3872  REM 

3874  REM  Check  that  initial  conditions  are  satisfied;  if  not,  then  start  over 
with  new  estimates  of  Tneutral  and  Tblood 
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Fblood  “"Ftfcl)  t  INPUT  "Enter  1  to  re 


3880  IF  INIT  )  0  THEN  PRINT  "Tr  «"T(11) 
start,  else  0"|  INIT 
3890  IF  INIT  >  0  THEN  GOTO  £130 
3900  REN 

3910  REM  Calculate  new  temperatures 
3950  REM 

3930  TIM  »  TIM  +  DT 

3940  FOR  I  «  1  TO  6  :  FOR  N  »  t  TO  3 
3950  T ( 10*I+N-10)  =  T(i0*I+M-10)  +  F <10*1 +N— 1 0 ) *DT 

3955  IF  T  ( 10* I +N— 1 0 >  <  TWAT  THEN  ORINT  "TEMP"  10*I-*-N-10"«"T  <  10*I+N-10>  s  END 

3960  NEXT  N 

3964  REM  Force  skin  tenioerature  to  aoaroacn  its  assigned  steady-state  value 

exponent lal ly;  when  Tskin  is  within  0. 0050  of  this  value,  zero  heat 
storage  assumed 

3965  TDIFF  =  T<  10*1-6)  -  BSTSK  :  IF  TDIFF  <  5.  000001 E-03  “HEN  TD  ( 10*1-6)  --  0(10* 

1-6)  -  E< 10*1-6)  -  BC< 10*1-6)  +  TD( 10*1-7)  :  GOTO  3967 

3966  T< 10*1-6)  -  SSTSK  +  TDIFF*EXP <F  < 10*1—6) *DT /TDIFF) 

3967  NEXT  I 

3970  T (61)  =  T  < 6 1 )  +  F(61)*DT 
3980  IF  TPRT  >  TIM  THEN  GOTO  4300 
3965  REM 

3990  REM  Print  results 
4000  REM 

4010  REM  Calculate  cardiac  output  (l/min)  and  bodv  ternoerature  (0) 

4020  CO  =  0  :  TB  =  0 

4030  FOR  I  =  1  TO  6  :  FOR  N  =  1  TO  4 

4040  CO  =  CO  BF ( l 0*1 +N- 10) /60 

4045  TB  =  TB  +  T < 10*I+N-10) *C < 10*1+N-10) /CB 

4050  NEXT  N  s  NEXT  I 

4070  REM  Calculate  total  evaporative  heat  loss  (kcal/h) 

4125  EV  =  E(l)  +  E (4)  +  E  <  1 1  ) 

4130  REM  Calculate  skin  ternoerature  and  skin  heat  flow  <kcal/h) 

4135  REM  Note  that  when  calculating  TDskin,  transition  occurs  when  Tskin  is 
within  0. 005C  of  its  steady-state  value 
4140  TSK  =  0  :  HFSK  =  0 
4150  FOR  1=1  TO  6 

416.0  TSK  =  TSK  +  T  ( 10*1-6)  *C  (  10*1-6) /CS 
4170  HFSK  =  HFSK  +  TD  ( 10*1-6) 

4180  NEXT  I 

4185  REM  Note  conversion  from  kcal/h  to  W 

4188  IF  1. 163*HFSK/SA  >  999  THEN  HFSK  »  999*SA/1.163 

4190  LPRINT  USING  "###.##  " ;  60*TIM, TB, 1. 163*MR/SA, 1. 163*EV/SA, T <61 ) , T ( 1 1 ) , T ( 1 ) 

, TSK, 1. 163*HFSK/Sft 

419S  IF  CODE  =  0  THEN  GOTO  4£00 

4195  FOR  I  =  1  TO  6  s  LPRINT  USING  "###.##  " ;  T < 10*1-9) , T ( 10*1-8) , T < 10*1-7) , T ( 1 

0*1-6)  :  NEXT  I 

4S00  TPRT  =  TPRT  +  DPRT 

4300  IF  TMAX  <  TPRT  THEN  END  ELSE  GOTO  £7B0 
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RESTING  WHOLE  BODY  COOLING 


HT  =  175.7  WT  =  79.  £  SG  =  1 . 062795  FBF  =  .1762  SA  -  1.9533GB 
QC  =  64. 9995  QM  =  8. 930359  QSF  -  5. 337824 
WC  =  £8.39902  WM  =  33.00977  WSF  =  17.79112 

Cot»e  Muscle  Fat  Skiri 


Weinnt  SWT  <kq> 


3.  485 

0.  384 

0.  465 

0.  £76 

14. 737 

18. 490 

8.  838 

1. 390 

2 .  2I8 

3.  476 

1. 210 

0.  498 

0.  269 

0.  078 

0.  186 

0.  192 

7.  146 

10. 504 

2.  977 

1 .  £36 

0.  445 

0.  078 

0.  £79 

0.  £45 

Rad i us 

R  ( cr«i ) 

9.  41 

9.  74 

10.  11 

10.  32 

8.  84 

13.  28 

14.  94 

15.  18 

2.  57 

4.  06 

4.  46 

4.  62 

0.  94 

1. 07 

1. 33 

1. 55 

3.  77 

5.  93 

G.  41 

6.  59 

1. 06 

1.15 

1 . 43 

1 . 63 

Therma 1 

Caoacitance  C 

<  kcal /C) 

2.  479 

0.  345 

0.  279 

0.  248 

10. 000 

16. 641 

3.  303 

1. 251 

1. 466 

3.  129 

0.  726 

0.  449 

0.  147 

0.  070 

0.  112 

0.  173 

4.  365 

9.  454 

1. 786 

1.  1 12 

0.  247 

0.  070 

0.  167 

0.  £21 

Basal  Heat  Production  QB  (kcal, 

1 3 . 656 

0.  104 

0.  139 

0.  083 

48. 503 

5.  002 

2.  652 

0.  417 

0.  647 

0.  940 

0.  363 

0.  150 

0.  075 

0.  021 

0.  056 

0.  058 

1. 994 

2.  842 

0.  893 

0.  371 

0.  124 

0.  021 

0.  084 

0.  074 

Basal  Blood  Flow 

BFB  (1/h) 

45.  00 

0.  12 

0.  17 

1.  47 

210. 00 

6.  00 

3.  18 

2.  17 

0.  78 

1.  13 

0.  44 

0.  52 

0.  09 

0.  03 

0.  07 

2.  01 

2.  39 

3.  41 

1. 07 

2.  94 

0.  15 

0.  03 

0.  10 

3.  07 
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Thermal  Conductance  TC  (kcal/C/h) 


1. 48 

6.  47 

8.  50 

0.  00 

2.  lc' 

4.  06 

10.  24 

0.  00 

3.  71 

7.  8) 

20.  87 

0.  00 

4.  85 

6.  1 6 

6.  89 

0.  00 

5 .  vi3 

1 1. 99 

35.  85 

0.  00 

6.  91 

8.  88 

9.  44 

0.  00 

Sepment  Surface  Area  A  <m*m) 

01  =  .133937  AS  =  . 5724  A3  =  .324942  A4  «  9. 346019E-02 
A5  «  .6630005  A 6  =  .1282577  AB  =  1.91599B 


T  neutra 

1  — 

air 

--■=  £< 

B.  85 

C 

Tbl 

ood 

=2  37  ■  S3 

C 

Twat, 

;er  —  £0  C 

' 

C0  =  24 

.  4 

3 

Cl  < 

=  65 

c 

C  “ 

5 

n 

L-l 

li 

• 

ui 

C4  = 

1.5  C5 

=  1 

C6  « 

0 

~i 

T 

c. 

T 

T  4 

37.  42 

35 

.  79 

35 

.  38 

i — *  • 

01 

37.  48 

37 

.  10 

35 

.  47 

34. 

07 

35.  15 

34 

.  58 

33 

.  85 

33. 

49 

35.  34 

35 

.  29 

35 

.  24 

35. 

17 

^6  •  42 

35 

.  70 

34 

.  76 

34. 

35 

35.  43 

35 

.  38 

35 

.  33 

35. 

£6 

TIME 

i 

B 

MR 

EV 

TBL 

TCR 

THD 

TSK 

HFSK 

0.  00 

34. 

97 

47. 

19 

5- 

83 

37.  25 

37. 

48 

37.  42 

34.  £8 

999. 

00 

6.  00 

33. 

42 

95. 

44 

9. 

04 

37.  £2 

37. 

45 

37.  33 

21. 03 

276. 

28 

12.  00 

33. 

02 

105. 

34 

9. 

70 

37.  17 

37. 

39 

37.  18 

£1. 00 

220. 

06 

18.  00 

32. 

73 

110. 

91 

10. 

06 

37.  1 1 

37. 

34 

37.  09 

21. 00 

198. 

72 

24.  00 

32. 

49 

1  15. 

36 

10. 

34 

37.  04 

37. 

27 

37.  03 

21. 00 

186. 

16 

30.  00 

32. 

£9 

120. 

64 

10. 

72 

36.  94 

37. 

19 

36.  95 

21. 00 

177. 

91 

36.  00 

32. 

13 

126. 

75 

11. 

10 

j6. 

37. 

08 

36.  85 

21. 00 

172. 

79 

42.  00 

31. 

99 

133. 

59 

11. 

53 

36.  71 

36. 

96 

36.  75 

21. 00 

169. 

65 

48.  00 

31. 

88 

141. 

00 

12. 

02 

36.  59 

36. 

84 

36.  64 

21. 00 

167. 

59 

54.  00 

31. 

79 

147. 

74 

12. 

52 

36.  48 

36. 

72 

36.54 

21. 00 

166. 

43 

60.  00 

31. 

71 

154. 

09 

12. 
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36.  38 

36. 
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